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This paper deals with the Autler - Townes doublet structure. Applied driv- 
ing and probing laser fields can have arbitrary intensities. The explanation is 
given of the broadening of doublet components with the growth of probing field 
intensity, which was observed in experiment. The effects of Doppler averaging 
are discussed. 
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I. INTRODUCTION 

The spectroscopy of three-level systems occupies a highly important place in nonlinear 
spectroscopy. In addition to saturation or power broadening of resonance , new remarkable 
effects appear: field splitting || and truly coherent processes (such as Raman scattering 
or two-photon processes). The most clear issue is the absorption spectrum of three-level 
systems interacting with two monochromatic waves within perturbation theory in intensity 
of one of the waves — probe field spectrum 0-^]. Some efforts were made to construct the 
theory where both waves are strong In ||, the coherences in the equations for the 

density matrix were excluded, and equations containing only populations were derived. As 
a result, the probabilities of transitions induced by fields were renormalized and expressed 
in terms of populations only. Unfortunately, our intuition in predicting the behavior of such 
systems is poor yet, it is difficult to foresee the result without solving the equations. The 
computation of three-level system with a great number of parameters (relaxation constants 
of levels, waves detunings and strengths) |§ is hard to analyze. 

However, power of light has increased since early experiments, as one needs to get better 



light conversion or more controllable setup. In experiment QTO | A-scheme interacting with 
two waves exhibits, at first sight, strange behavior. When the intensity of comparatively 
weak wave increases, the components of Autler - Townes doublet broaden in its absorption 
spectrum; finally the doublet is transformed to a single line. 

A three-level system with two strong fields can be a part of the scheme of four-wave 
resonant mixing. The study of this area is progressing rapidly, since it gives a hope to 
obtain coherent CW short-wave radiation. To understand the whole picture it is useful at 
first to consider the part of a system that interacts with strong fields. 

The aim of the present paper is to analyze qualitatively the effects arising from increasing 
the power of weak wave in A-scheme. In what follows we discuss the general features of 
multi-level systems and the possibility of reducing them to simple two-level systems. We 



use the probability amplitudes approach and the "dressed" state concept (see JTTIJT2 ] ) to 
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understand the spectra. The energies of "dressed" states are identified with the poles of the 
Green function of master equations in cu-representation. We are not interested in intrinsic 
relaxation processes in an atomic system and this allows us to avoid introducing the density 
matrix. It reduces the system of N 2 equations for the elements of the density matrix to the 
system for N probability amplitudes (here N — 3). 
The paper is organized as follows: 

In sections P|,pi| the equations of motion and their Green function are described. In 
section |V| the energies of "dressed" states are discussed. Section |V| is devoted to the calcu- 
lation of probe field spectrum in a four-level system — linear response of a three-level one; 
it shows the benefits of Green function. Section [VT| presents the absorption spectra of the 
strong field. In section [VI]] we briefly consider the question of the Doppler broadening of 
spectra. In section |V111| we finally summarize essential aspects of the paper. 



II. BASIC EQUATIONS 

Consider the atomic system consisting of three excited levels and its resonant interaction 
with two electromagnetic waves with frequencies u>i, u 2 (see figure 1). Hereafter we assume 
that the wave amplitudes are not changed essentially by interaction (optically thin media). 
Thus, we do not consider here the effects of light propagation. One can think that there is a 
thin layer with absorbing atoms and light travels trough it almost without absorption. Here 
we also do not take into account the motion of atoms, i.e., we do not consider the effects 
connected with inhomogeneous broadening of absorption spectra. Then the time evolution 
of atomic wave function is described by the Schrodinger equation (h — 1) 




-G*e iwi * 

e 2 -±t 2 -G* 2 e^ 1 

de-*" 1 * -G 2 e-^* E 3 - ir 3 




(1) 
(2) 



Here 7{q is the atomic Hamiltonian, 7^ int describes the interaction with light; a iy E iy Tj 
are the probability amplitude, the energy, and the decay rate of state Gi is the matrix 
element (3| — 'H\nt\i)^ luJit of interaction. Equation (0) is written in resonant (or rotating 
wave) approximation, which is valid if Gi, Ui — E 3 + E^ -C E 3 — E iy % = 1,2. The frequency 
Qi = u>i — E 3 + Ei is the detuning of the wave i from resonance. 

Since there is no closed cycles made up from wave fields (see figure 1(a)) one can bring 
(0) to the form without explicit dependence on time: 



dt 



A i = a i ^ Ei - n **, i = l,2, A 3 
f A t \ /fii-iri 

o n 2 - ir 2 



a 3 e iE3 *, 



(3) 

A 2 I = I Q 2 - |r 2 -GI | A, | . (4) 

Moreover, fields Gi can be thought of as real. Nevertheless, sometimes we will write G* 
instead of Gi. It helps to determine the interaction process which corresponds to the written 
expression. 
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In spectroscopy, the observed quantities are the powers absorbed from the waves. Light 
field induces transitions between atomic levels, and absorption is connected with the change 
of levels population. The time derivative of population on level 2 is equal to 

A|A 2 | 2 = -r 2 |A 2 | 2 - iG 2 A 2 Al + iG* 2 A* 2 A 3 . (5) 

The first term in r.h.s. of (|5|) is due to the decay of the state |2). Other two terms are 
responsible for the absorption of the 2nd wave. 

The formulation of the problem of light absorption is the following: at the moment to 
the atom is excited to the state A = Ao- Then the time evolution of A is governed by (|]), 
and the power absorbed from the 2nd wave by this atom is equal to 

/•OO 

p(a , n 2 , t ) = uj 2 at (iG 2 a 2 a; - iG* 2 a;a 3 ) . (6) 

Because equation (|j) is time-independent, the power P does not depend on to. Then one can 
set to = 0. In addition, P is a quadratic functional of A(t), and equation (^) is linear. Thus, 
P can be written as P(Aq,Q 2 ) = uj 2 AqV(Q 2 )Ao, where ^ is the Hermitian conjugation. 
The matrix V contains complete information about the absorption of the 2nd wave, i.e., 
using V one can calculate P with any initial excitation A . For example, V 22 is equal to the 
probability of absorption of a photon from the 2nd wave if the atom was initially excited to 
level 2. Because P is real, the matrix V is Hermitian. 



III. GREEN FUNCTION 

Let us consider equation (|) with the initial condition Kit = 0) = A : 

i^A = HA + i8(t) A , A(t)=0 if t < 0. (7) 

Since this equation is time-independent, it is convenient to consider its Fourier transform, 
i.e, its energy representation (h = 1) 

A(t) = J ^ e-^A,,, ujA w = HA U + iA . (8) 

It is a set of linear algebraic equations. Its solution is equal to 

A w = i(u - Hy'Ao = W(u)A . (9) 

The matrix T>(u>) is the Green function of (£|) in cj-representation; it is also the resolvent 
of the Hamiltonian H . The matrix element T>ij(u) describes transitions i <— j induced by 
fields. As a function of complex variable u, the matrix T>(u) has poles at the eigenvalues 
\j, I = 1,2,3 of the Hamiltonian H, i.e., at the quasi-energies of the stationary states of 
"atom-l- field" system ("dressed" states). The energy A of a "dressed" state is found from 
the equation Det (A — H) = 0: 
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(a - ^ + |r x ) (a - n 2 + |r 2 ) (a + |r 3 



a - ^ + |rij |G 2 | 2 + (a - n 2 + \v 2 ) idi 2 . (10) 

Each A/, J = 1, 2, 3 has negative imaginary part Im Xj < corresponding to the decay rate 
of the "dressed" state |/)). In other words, P(o;) has poles only in the lower half-plane. 

If the fields G\, G 2 are weak, then there is a hope that the evolution of the system is little 
different from the evolution without fields. One can use this "bare" evolution as a starting 
point and expand the Green function T>(Gi,G 2 ) into a power series in the neighbourhood 
of Gi = 0, G 2 = 0. We denote 
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(11) 
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Then H = E + G. One can write the following expansion: 

V = D + DGD + DGDGD + ... = D + DGV, (12) 

where the argument of T> and D is equal to lu, and D{uj) = (u — E)~ x . This series is an 
analog of (1 — a;)^ 1 = 1 + x + x 2 + . . .. "Bare" Green function D(u) has poles at the energies 
of "bare" states Qi — ^r 1; Q 2 — ±T 2 , — |r 3 . The summation (p!2|) shifts the poles. 

The diagonal element T>a is the Green function of atom in the state It also has poles 
at Xj. There is an admixture of all dressed states |/)) in In the limit of weak fields T>a 
has obvious pole. If the fields induce transitions between the states |i) and \j), then the 
atom in the state \i) lives in the state \j) for some time. Then the pole corresponding to \ j) 
appears in its Green function. When the fields become stronger these poles shift. 

Now we can write a simple expression for the matrix V(Q 2 ) describing the absorption of 
the 2nd wave. 



P(A ,n 2 )=u 2 dtA\t)PA(t), iP 
Jo 



/0 

G* 2 | . (13) 
V0 -G 2 



Vector A(t) can be rewritten in terms of the Green function. We obtain 



V\uj)PV(uj). (14) 



&uji&uj 2 V^(uji)PV(tj0 2 ) f did ^ 



(2vr) 2 lo 2 - 00! - iO J 2vr 

The correction iO in (|14|) is due to causality. The physical determination of the integral over 
t consists in multiplying by a factor e~ £t with e — > +0. Since to (as well as fli and Q 2 ) is real 
on the contour of integration, in the integral we can replace T>\uj) by T) c (oj). The operation 

■ c is defined as follows: - c = f -\ but during this operation the frequencies u>, Qi, i — 1, 2 are 
thought to be real, i.e., they need not be complex conjugated. The operation - c changes only 
the sign of Ts. Thereafter the integrand in (|T4]) is an analytic function of uj, i.e., it does 
not depend on u*, Q*. Moreover, it is rational; one can apply residue theory to calculate 
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the integral. In its turn, V{Vt 2 ) is an analytic function of Q 2 and the poles of V correspond 
to the resonances in the 2nd wave absorption spectrum. The real part of the pole position 
gives the frequency of the resonance, the imaginary part represents the resonance width. 

The matrix H has right and left eigenvectors: H\I)) = Aj|/)), ((I\H = ((/|Aj. They can 
be normalized so that ((I\J)) = Sjj. Then 

H=f:\i))Xi({i\, v(u)=j:\i}}^—((i\. (is) 

1=1 1=1 w - A/ 

If V(Q 2 ) has a pole at fl 2 = fi, then the integral (JTJ]) for Q 2 — ^ is divergent. The contour of 
integration bypasses the poles of £>(u;) and T> c (uj) in different directions. Then the integral 
([14]) may be divergent only if T> and T> c have at least one common pole Xj = Xj. When Xj 
and Xj merge together they squeeze the integration contour, and the integral ([H]) acquires 
singularity. Time evolution of the dressed state |/)) is described by the exponent e~ lA/ *. The 
case where Xi is close to A j corresponds to the resonance between the waves and the dressed 
states |J)), | J)). The absorption matrix V(Q 2 ) contains resonance denominators of the form 

A/ " A';. 



IV. DRESSED STATES 

Here, for simplicity, we neglect relaxation constants. In fact, the resonance widths under 
power broadening depend on the ratio of Ts even when the latter are small. Therefore, the 
contents of this paragraph is useful for the case of equal relaxation constants (or, maybe, 
when they are comparable to each other) or for finding resonance positions in the absorption 
spectrum of auxiliary probe field. The energies of the dressed states satisfy algebraic equation 
(PH|). Let us write it in the form 

[A(A - Qi) - |G X | 2 ](A - n 2 ) = (A - fii)|G 2 | 2 . (16) 

If we consider the energy levels of the system "atom+field" , then we will see the triples of 
close levels. For weak fields these levels are: \3)\rii,n 2 ) (atomic state |3), n x quanta in the 
1st field, and n 2 quanta in the 2nd one), |l)|n! + l,n 2 ), and \2)\rii,n 2 + 1). If we consider 
the system within perturbation theory, we use the atomic basis of states. They are close 
in energy; in resonant approximation the field-induced transitions occur only inside each 
triple. The "dressed" energy levels are found from the secular equation for each triple (|10|) . 
In resonant approximation the triples are independent. 

The resonances in the spectra occur when two eigenvalues A are close to each other. The 
minimal distance between them gives the order of resonance width. The general idea is 
that one can consider two close As separately from others and write the quadratic secular 
equation for them. These two eigenvalues are pushed apart by the intermixing — the general 
scenario due to the Hermitity of interaction. 

When G 2 < Gi, ^i,2> A is found from the condition of vanishing of l.h.s. in (p"6|): 



A 1)2 = (Hi ± (Of + 4|C 1 | 2 ) 1/2 )/2, X 3 = n 2 . (17) 

First two dressed states result from |1) and |3) by the 1st wave field splitting. The third 
one corresponds to |2). The 1st wave makes the distance between Ai and X 2 not less than 
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2|Gi|. Then the resonances in the 1st wave absorption spectrum undergo power broadening 
— their width becomes proportional to |G*i|. 

When O2, i.e., A3, is close to Ai, we can consider Ai and A3 separately from A2. Denoting 
^2 = Ai + Q, A = Ai + A, we obtain 

Mi,2 = U 1 * ,^ . ^ , 2 ,i/ 2 V A(A-fi)^M 1 |G 2 | 2 . (18) 



(fi? + 4|Gi 



The energies Ai and A3 are split by the 2nd wave, whose intensity is multiplied by the factor 
Mi. The case A3 ~ A 2 is handled in a similar way, one should only substitute A 2 and M 2 
for Ai and M±. The coefficients M 12 were called "memory factors" or "correlation factors" 
|5|Jl3[|. One can diagonalize the Hamiltonian H neglecting the 2nd field G 2 and turn from 



bare states |1) and |3) to the states dressed by the 1st field having energies Ai >2 . Then 



Mi i2 |G 2 | are simply the matrix elements of interaction between these dressed states and 
bare state |2), induced by the 2nd field. It is well known that while the 2nd wave is weak 
(probe), its absorption spectrum P(Q 2 ) has two resonance lines (Autler - Townes doublet) 
||. In this case their width is determined by the relaxation constants. From the behavior 
of the dressed states we conclude that these two lines will undergo power broadening by the 
2nd field if it becomes stronger. 



V. PROBE FIELD SPECTRUM 

We consider the scheme illustrated in figure 1(b) in the case where only the 4th level 
is excited (i.e., A = (0,0,0,1)). We assume the fields G M and G u to be weak. Thus, 
only the 1st order of perturbation theory in G^ v is needed. In the main order, A4 = 
exp (— r 4 t/2 + i(f2i + flfj,)t) ■ 9(t), where At = 0,4 exp (i(i^ 4 — Qi — fl^t). The probability 
amplitudes Ai, i = 1,2,3 have the form 

A(t) = J dt (G^Vait - t ) + G u e- i£t V l2 (t - t j) A 4 (t ) 

du iwt ( GyVjiju) G u V l2 (u) \ 

2ir \F 4 /2 + + ty, - u) T 4 /2 + i(Q 2 + Q„ - u) J ' 1 ' 

where e = Q 2 + Q u — Qi — and T>{t) = J du> e~ llJJt T>(uj) /2tt . The power absorbed from the 
field G^ is equal to 



/•oo r 

P M = 2uv Re iG; j dt A\A X = 2u„ Re iG^ J 



du 1 



2vr T 4 /2 - i(fii + Vt^ - u) 



+ 



^ r 4 /2 + i(fi x + ^ - 00) r 4 /2 + i(n 2 + n 



UJ 



^n eiG; (^^ + 9pm, (20) 

where u;* = fl± + + ir 4 /2. The expression for the absorption P At (fi M ) linearly depends on 
the Green functions T>\\ and T>i 2 , since we considered the linear response of the three- level 
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system. Note that if the initial excitation goes to levels 1, 2, 3, then the result will contain 
the product of two Green functions. If one tries to calculate the power absorbed from the 
field ( pop using ([14]), where two Green functions and the matrix P are taken for the 
four-level system, then one of these Green functions (T> or T) c ) will be equal to P44 ~ -D44. 

The second term in brackets in (^) corresponds to parametric process. There is a 
closed cycle made up from fields /i-l-2-z/ with the whole detuning e. This contribution to 
the spectrum P^e) has the width T4. However, this is only the gain in the thin media 
approximation (or the parametric instability increment). If the parametric process leads 
to instability, then the emission spectrum of thick media has different width, which can be 
much smaller. 



VI. SPECTRUM V(n 2 ) 

Here we consider the case of small equal relaxation constants I\ = T 2 = T 3 — » 0. We 
focus only on the diagonal elements of the matrix V, which give the absorption power when 
the initial excitation goes to the real atomic levels. For real Gj, j = 1,2, the components 
of the eigenvectors of H are also real. From ([T4]) and flT5| ) immediately follows 

w ^ E fMI . (21) 

i±j X J ~ X i 

As a function of As, ([H]) is a rational expression symmetric with respect to permutations. 
Then it can be easily expressed in terms of the coefficients of i.e., waves detunings and 
intensities. Note that TrP is equal to zero, since ((J\i) = &u (in the case of equal 

populations of levels the interaction is absent). The power absorbed from the 2nd wave 
when the excitation goes to level 2 is given by 

V 22 = 2G 2 2 P/Q, h = Gl, h = G\, I 12 = G\ + Gl (22) 

p = (4Ji 2 + nl)(i 12 + (fix - n 2 ) 2 ) - ml - j 2 [3/ 12 + m\ - ioo x q 2 + m 2 2 ], (23) 

q = (4/ 12 + nj)(i 12 + (q 1 - n 2 )n 2 ) 2 + / 2 (fi! - n 2 ) 

■[2(Qt - 2ft 2 )(9/ 12 + (fix + Q 2 )(2Q 1 - Q 2 )) - 27I 2 (Q! - 2 )]. (24) 

Consider the simple case Q\ = 0. One can measure frequency in \Gi\ units, i.e., set 
\Gi\ = 1 and obtain 

V ^ GI + 5G 2 + A + ml 

V22iSh) ~ 2G2 A{Gl + 1)3 + (Gl + 20G| - 8pl + 4fif (25) 
The poles in the spectrum are situated at the points 

a ^ ± ^^^ ± M^f!lf, (26) 

(for the poles positions in other simple cases see ||). The asymptotics for small and large 
G 2 are the following: 
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_ / 1 ± V2iG 2 - Gl/ A, G 2 « 1, . s 

U2 *~ ± \2\G 2l iGl/2, G 2 »l. 1 /J 

The poles ^2* = ii^l/^ f° r G2 ^ 1 are inessential, since their contribution to the spectrum 
is negligible. The positions of the poles Q 2 * = ±2iG 2 correspond to the usual power broad- 
ening by the 2nd wave or to saturation HJ. The resonance positions Q 2 * ~ ±1 for G 2 <C 1 
are the positions of Autler - Townes doublet components 0. The addition v^2iG2 is the 
power broadening of these components by the 2nd field. The term — G 2 /4 is the nonlinear 
shift of the resonance produced by this field (see |A]). 

When G\ — 8 the poles in the spectrum P(VL) merge together: fi 2 * = ±3 3 / 2 i, V 22 = 
4(104 + 3^2)/(27 + Ql) 2 (see figure 4 (a)). At G 2 = V8 all dressed state energies are equal 
to each other: Ai )2 ,3 = ^2*/3 = ±v3i. This situation seems to be destroyed when fii ^ 
(see figure 4(b)) or relaxation constants are not negligible. 

In figure 3 one can see that the power broadening of Autler - Townes doublet components 
is proportional to the amplitude of the 2nd wave. When Qi 7^ the broadening is asymmetric 
because the rates of interaction between level 2 and the states dressed by the 1st field are 
not equal to each other (Mi 7^ M 2 ). When the 2nd field becomes stronger the resonances 
come closer to each other. When G 2 G\, Q± the poles in the spectrum are situated at the 
points 

ft 2 * = ±2iG 2 , ft 2 * = 7T~~TJ~<v ' ( 28 ) 

The first pair of poles corresponds to the simple resonance between levels 2 and 3 taking 
into account the saturation caused by the 2nd field. If we decrease the intensity of the 
2nd wave, then this resonance will be transformed to the component of Autler - Townes 
doublet Q 2 = (fij — (f2 2 + 4|Gi| 2 ) 1 / 2 )/2. The second pair of poles corresponds to the other 
component of the doublet (if the 1st field is weak, then this line is the resonance between 
the virtual level and level 2, i.e., two-photon resonance). The position of this resonance for 
G 2 ^> Gi,Qi (the real part of the pole position) has the sign opposite to that of f^. Note 
that in figure 3(b) the narrower part of the spectrum goes to the left. For extremely high 
values of G 2 the second part of the spectrum becomes wider than the first one (the width 
becomes proportional to Gl 1 , see dots in figure 4(b)). 

The fact that the width of one of the two resonances is proportional to G\ for G 2 3> G\, Qi 
corresponds to the following: the 2nd field has the amplitude G 2 , which is much greater than 
the distance (f2 2 + 4 1 Ci | 2 ) x / 2 between the states dressed by the 1st field. Then the 2nd field 
feels these two states as only one state; it cannot see the energy structures with resolution 
higher than \G 2 \. As the system becomes similar to an ordinary two- level system, only one 
resonance with the power width 2\G 2 \ should remain. There are different ways to kill the 
2nd resonance: to increase its width or to decrease its amplitude, both ways being used here. 
Mathematically, the fourth-order algebraic equation for f2 2 * is transformed to the second- 
order algebraic equation, i.e., the coefficients of f^* and f^* are comparatively small. It is 
well known that the algebraic equation with small leading coefficient e has at least one root 
which is large in parameter 1/e. 
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VII. DOPPLER BROADENING OF SPECTRA 



One can take into account the Doppler broadening of absorption spectra by substituting 
Qi — kiv and Q 2 — for Qi and Q 2 and integrating the result (|22l) with some distribution 
of particles p(v). Consider the case of copropagating waves, where only a longitudinal 
projection vu of velocity is needed. When the Doppler width is infinite, the answer for the 
spectrum is the integral of the rational function p±(y \\)/p§{v\\) over velocity, where p n (x) is 
the nth degree polynomial. Some simplification takes place when ki = k 2 or ki = for 
some i = 1, 2. In the general case the expressions for velocity-averaged ^(^2) are difficult to 
derive and analyze. We will discuss the effects of the Doppler broadening only qualitatively. 

When the field G 2 is weak and k 2 < ki, the velocity-averaged spectrum V(Q 2 ) has two 
narrow asymmetric resonance lines coming from turning points |14| of frequency branches 



Q 2 ,(v ]{ ) ~ M|| + (fix - hv\\ ± ((fii - k lVll ) 2 + 4|Gi| 2 ) 1/2 )/2. (29) 

The turning point is the extreme point of velocity dependence of resonance frequencies fl 2 *, 
i.e., at the turning point v tp we have (df2 2 */df ||)L|=i> t p = 0- The spectra of this type are 
well known (see, e.g., ]15|]). For the universal shape of asymmetric line due to the turning 
point see ||16|| . When G 2 <C T, the characteristic width of this resonance line is of order T. 
When G 2 S> T, the resonance line width is determined by the power broadening by the 2nd 
wave. Nevertheless, while G 2 <C G± one can use the expression (|29|) (where the influence of 
G 2 on the resonance position is neglected) for f2 2 *. The rough feature of the spectrum - 
two narrow asymmetric resonance lines — remains, but the width of these lines depends on 
G 2 (see figure 5). 

When Qi and G\ are compared to the Doppler width, an "isolated peak" appears in the 
spectrum. Its width depends on fix, G\ (see |T^fi7[1 ) and is larger than G 2 in the general 
case. The dependence of the width on Q 1 and Gi due to inhomogeneous broadening is not 
a subject of the present paper, but at certain conditions the width decreases and becomes a 
value of order G 2 . 

When G 2 ~ G\ ^> T, the absorption spectrum does not contain narrow lines, i.e., all 
resonance widths are of order or \G 2 \. 

If k 2 > ki, then the frequency branches (|29|) have no turning points, and the absorption 
spectrum V(Q 2 ) has no narrow resonance lines even when the 2nd wave is weak. If we 
increase the intensity G\ the narrow lines will not appear. 



VIII. CONCLUSIONS 

Let us summarize the description of the three-level system resonantly interacting with 
two strong monochromatic waves. 

The expression for the spectrum contains resonant denominators Xj — A}, where A/ is the 
energy of 7th "dressed" state. For iV-level atomic system the dressed state energies satisfy 
the algebraic equation of iVth order which is the secular equation for N levels of the system 
"atom+field" that are close in energy. When two dressed states are close to each other and 
all other states are far from them in energy, the secular equation can be reduced to the 
quadratic secular equation for these two states. The part of the atomic system consisting of 
these two dressed states behaves as an ordinary two-level system. 
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If some connected component of strong fields covers the whole atomic system with N 
states, then all N "dressed" states will be far from each other. As a consequence, all 
resonances in the spectra will be wide. Atomic system can be covered by strong fields, 
but strong fields can be disconnected (e.g., fields 1 and v in figure 1(b) are strong). Then 
the matrix elements between the states dressed by field 1 and by field v will be small and 
resonances will be narrow. 

When the 2nd wave is weak we can think that it couples level 2 (see figure 1) and levels 
1, 3 split by the 1st wave. The distance between the split levels is much greater than the 
amplitude of the 2nd wave; these two transitions cannot be resonant simultaneously. Each 
transition can be treated as a two-level system. The effects of power broadening by the 2nd 
wave appear. When the 2nd wave is very strong, only the two-level system coupled by it 
remains and Autler - Townes doublet becomes blurred. 

Indeed, the power width in a two-level system depends on the ratio of relaxation constants 
T of levels even when the latter are small. Thus, the power width of Autler - Townes 
component, due to the intensity of the 2nd wave, should also depend on the ratio of Ts. The 
expressions will be similar to the well-known ones for a two-level system: one should take 
T 2 and the decay rate of the state "dressed" by the 1st field. 

In experiment |]TD| the 1st wave detuning Qi was equal to zero and two components had 



the same power width. It is interesting to observe the asymmetric broadening when ^ 0. 

When the 2nd wave is generated in the media, self-consistent problem taking into account 
energy losses should be solved. The doublet gives two frequencies of generation, which merge 
together at a certain rate of reflection losses. 

The author is grateful to E V Podivilov and S A Babin for useful discussions. The present 
paper was partially supported by Soros Foundation (gr. a98-674), by INTAS (gr. 96-0457) 
within the program of ICFPM, and by RFBR (gr. 96-02-19052). 



APPENDIX A: RESONANCES IN PROBE FIELD SPECTRUM 



Consider the n-level system with dressed energies Ai, A 2 , . . . , \ n and probe field spectrum, 
where the probe field G with detuning Q resonantly interacts with level |0) and the n- level 
system. The dressed energies A of the whole system can be found from the Hamiltonian: 



H 



X-Q 



f M[G 
rW x G 



r MlG JW 2 G 
Ai 
A 2 



\VM~G 

Mi 





M 2 



f M~ n G\ 





A„ / 



The resonances are situated at ft, 
|) has multiple root if 

= Ai ± 2r/MiG - 



A — Ai A — A 2 

Ai, i 



+ . . . + 



G 



(Al) 



(A2) 



1,2, 



A -A, 

n. When Q is close to Ai the equation 



M 2 
Ai-A : 



+ . . . + 



M n 
Ai — A T 



G 2 + 0(G 3 ). 



(A3) 
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The Green function Vqo(uj) has the form 

V~\u) = (( W - H)^)- 1 = a; - O - £ (A4) 

■ U! — Aj 

The mass operator oj — Q — Vqq is the sum of loops through levels i = 1,2, ... ,n. When 
Q ~ Ai the loop i — 1 is much greater than others. This loop gives the power broadening 
of the resonance Q ~ Ai. Other loops renormalize the detuning, which gives the shift of the 
resonance. Note that the shift is of order G 2 / (Ai — Aj), which is much smaller than the width 
of the resonance ~ G. The shift of the resonance in (|A3|) can be important only when the 
field G is not very weak. Then it gives qualitatively correct answer. When the probe field 
G is weak enough, the shift can be neglected and the resonances Q ~ A» live independently. 



APPENDIX B: RESONANCE POSITIONS IN SPECTRA 



Although it is not very difficult to calculate the integral (|T4]), some properties of V{Vt 2 ) 
can be elucidated by simple algebraic manipulations. It should be noticed that the same 
result (also without integration) can be obtained using the density matrix formalism. 

The resonance position f2 2 = f2 in the 2nd wave absorption spectrum can be found 
from the condition that one pole of T> coincides with at least one pole of T> c . Consider the 
following combination of eigenvalues: 

Z = H(\j - A}) = Det^(£ c ) = -Detf 6c (H), (Bl) 
i, j 

where /^(A) = Det (A — A). The resonance occurs when Z = 0. As a function of As and 
A c s, Z is symmetric with respect to permutations, so it can be easily expressed in terms 
of detunings, fields, and relaxation constants. In the general case the expression for Z is 
bulky, hence we will consider some particular cases. We assume that I\ = T 2 = T 3 = T; 



n t = n 2 = n (iD, ^i = o O), n 1 = r = o (pp. 

Z = -T 3 Z, h=Gl h = Gl, I 12 = Gl + G 2 2 , (B2) 

z = (n 2 + 4/ 12 + r 2 ) (r 2 n 2 + (i 12 + r 2 ) 2 ) , (B3) 

z = (4/ 12 + r 2 )((^ 2 - / 12 - r 2 ) 2 + 4T 2 n 2 2 ) - i 2 n 2 2 (4n 2 2 - 36/x - 9/ 2 ), (B4) 

Z = 4/ 12 (fi 2 - J 12 ) 2 - J 2 fi 2 (4fi 2 - 36/i - 9I 2 ). (B5) 
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FIGURE CAPTIONS 

FIG. 1. 3-Level system interacting with two electromagnetic waves (a) and its testing 
by two weak (probe) waves (b). 

FIG. 2. Dressed state energies A(f2 2 ), G x = 1. One can see the asymptotics Ai )2 , A 3 (jPTp. 
When A3 ~ Ai,2 the energies A repel each other with the rate (iWi^) 1 / 2 !^! (see ([18])). Figure 
2(a): fii = 0,' G\ = 1/32 (solid curve), G\ = 1/2 (dashed curve). Figure 2(b): fi a = 4, 
G\ = 1/8. When tt x ^ we have M x ^ M 2 , the splitting by the 2nd wave is assymmetric. 

FIG. 3. Absorption power P 22 (fi 2 ,G 2 ) at G x = 1, Vt x = (a) and Vt x = 4 (b). The 
amplitude of the 2nd field G2 varies from 0.2 to 1.4 with the step 0.2 (a), and from 0.5 to 3 
with the step 0.5 (b). 

FIG. 4. Pole positions f2 2 * = (Ref2 2 *, Imf2 2 *) for different values of G2, G x — 1, Q x — 
(a) and fl x — 4 (b). The dots in figure 4(a) correspond to G2 = \fn, n — 0, 1, . . . , 8; in 4(b) 
— to G 2 = 0, 1, . . . ,5. 

FIG. 5. Velocity- averaged spectrum P 22 (fi 2 ) (numerical calculations), G x — 1, Q x — 0, 
k x = 1, = 0.8, G2 varies from 0.2 to 1.4 with the step 0.2. Uniform velocity distribution 
or infinite Doppler width. 
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